Modelling of groundwater flow in heterogeneous porous media by finite element method     

Hassan Smaoui  Lahcen Zouhri  Abdellatif Ouahsine  Erick Carlier 《水文研究》2012,26(4):558-569

The HySuf‐FEM code (Hydrodynamic of Subsurface Flow by Finite Element Method) is proposed in this article in order to estimate the spatial variability of the transmissivity values of the Berrechid aquifer (Morocco). The calibration of the model is based on the hydraulic head, hydraulic conductivity and recharge. Three numerical tests are used to validate the model and verify its convergence. The first test case consists in using the steady analytical solution of the Poisson equation. In the second, the model has been compared with the hydrogeological system which is characterized by an unconfined monolayer (isotropic layer) and computed by using PMWIN‐MODFLOW software. The third test case is based on the comparison between the results of HySuf‐FEM and the multiple cell balance method in the aquifer system with natural boundaries case. Good agreement between the Hydrodynamic of Subsurface Flow, the numerical tests and the spatial distribution of the thickening of the hydrogeological system is deduced from the analysis and the interpretations of hydrogeological wells. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

Simulation and Validation of the Aerosol Optical Thicknessover China in 2006     

ZHANG Hu  ZHANG Min  CUI Zhenlei  WANG Yuesi  XIN Jinyuan 《Acta Meteorologica Sinica》2012,26(3):330-344

  相似文献   

Crustal complexity in the Lachlan Orogen revealed from teleseismic receiver functions     

F. R. FONTAINE  H. TKALCIC  B. L. N. KENNETT 《Australian Journal of Earth Sciences》2013,60(3):413-430

There is an ongoing debate about the tectonic evolution of southeast Australia, particularly about the causes and nature of its accretion to a much older Precambrian core to the west. Seismic imaging of the crust can provide useful clues to address this issue. Seismic tomography imaging is a powerful tool often employed to map elastic properties of the Earth's lithosphere, but in most cases does not constrain well the depth of discontinuities such as the Mohorovi?i? (Moho). In this study, an alternative imaging technique known as receiver function (RF) has been employed for seismic stations near Canberra in the Lachlan Orogen to investigate: (i) the shear-wave-velocity profile in the crust and uppermost mantle, (ii) variations in the Moho depth beneath the Lachlan Orogen, and (iii) the nature of the transition between the crust and mantle. A number of styles of RF analyses were conducted: H-K stacking to obtain the best compressional–shear velocity (V _P /V _S) ratio and crustal thickness; nonlinear inversion for the shear-wave-velocity structure and inversion of the observed variations in RFs with back-azimuth to investigate potential dipping of the crustal layers and anisotropy. The thick crust (up to 48 km) and the mostly intermediate nature of the crust?mantle transition in the Lachlan Orogen could be due to the presence of underplating at the base of the crust, and possibly to the existing thick piles of Ordovician mafic rocks present in the mid and lower crust. Results from numerical modelling of RFs at three seismic stations (CAN, CNB and YNG) suggest that the observed variations with back-azimuth could be related to a complex structure beneath these stations with the likelihood of both a dipping Moho and crustal anisotropy. Our analysis reveals crustal thickening to the west beneath CAN station which could be due to slab convergence. The crustal thickening may also be related to the broad Macquarie volcanic arc, which is rooted to the Moho. The crustal anisotropy may arise from a strong N–S structural trend in the eastern Lachlan Orogen and to the preferred crystallographic orientation of seismically anisotropic minerals in the lower and middle crust related to the paleo-Pacific plate convergence.  相似文献   

An approach to determining the thicknesses of shear bands with an echelon‐crack structure     

G. S. Wang  L. W. Kong  C. F. Wei 《国际地质力学数值与分析法杂志》2013,37(6):618-640

In this paper, the internal structure of shear band is investigated, and a model of the shear band with an echelon crack structure is developed. The model assumes the shear band to be composed of two conjugate sets of echelon cracks, such that the smaller echelon cracks are embedded in the space of the larger ones. The additional strain induced by the echelon cracks and the anisotropic development of the compliance tensor in the shear band zone are analyzed. The critical crack density at the onset of shear band is obtained by applying the bifurcation condition. Deviating from previous approaches, the new procedure evaluates the thickness of shear band based on the geometrical characteristics of echelon crack arrays and the failure probability of grain boundaries in the longitudinal direction at the onset of shear band. Parametric analysis shows that grain size, internal friction angle, dilation angle, and failure probability of grain boundaries are the dominant factors that account for the shear band thickness. The calculated results are consistent with the experimental data available in the literature. The model soundly explains that the measurements of the shear band thickness are generally scattered, ranging from 4 to 30 (or even more) times the grain size. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

1961-2020年郑州主城区大气自净能力变化特征分析     

刘美  姬兴杰  田力  丁亚磊  左璇 《气象与环境学报》2022,38(3):93-100

利用郑州市主城区1961—2020年气象观测资料和2014—2018年空气质量监测数据,分析了郑州主城区大气自净能力指数的长期变化趋势与影响因子以及2014—2018年主城区大气自净能力与PM_(2.5)的关系。结果表明:郑州主城区大气自净能力指数30 a气候均值为4.42 t·(d·km^(2))^(-1),春季大气自净能力最强,为5.20 t·(d·km^(2))^(-1);秋季大气自净能力最弱,为3.88 t·(d·km^(2))^(-1),不利于对大气污染物的清除。1961—2020年郑州主城区大气自净能力呈显著的减弱趋势,其中1969年最强为6.85 t·(d·km^(2))^(-1),2020年最弱为3.06 t·(d·km^(2))^(-1)。影响因子中,1961—1980年混合层厚度与大气自净能力指数呈正相关;日平均风速≥2.5 m·s^(-1)的日数和小风日数与大气自净能力分别呈正、负显著相关;大气自净能力指数与降水日数显著相关,2015年后偏强降水日数的增加对大气自净能力在同时期的增强有一定影响。此外,研究还表明主城区大气自净能力和PM_(2.5)浓度存在显著的负相关,说明大气自净能力强时,对应的PM_(2.5)浓度低,环境空气质量趋好。  相似文献   

谢一矿小型层滑构造特征对煤厚变化影响分析   总被引：1，自引：0，他引：1  

殷寿明  黄晖 《安徽地质》2005,15(3):190-192

通过分析层滑构造的基本形式、分布特征和在煤层中的展现形态,指出层滑构造对煤层厚度变化的影响特征,对煤厚的变化做到定性和定量预测.  相似文献   

大陆岩石圈强度研究     

熊熊  王英杰  李永东 《大地测量与地球动力学》2022,42(5):441-450

岩石圈强度或有效弹性厚度(Te)控制着岩石圈对长期负载的响应及其演化过程和空间构型,包含丰富的地球动力学信息,对解译地壳、岩石圈乃至地幔介质的力学性质和动力学过程及机理具有重要的科学意义。本文拟回顾从均衡概念的提出至今关于岩石圈强度研究的发展历程,重点从历史角度梳理岩石圈强度研究的沿革脉络、各种方法发展的逻辑关联,并对岩石圈强度与壳幔耦合、岩石圈有效弹性厚度与地震孕震层厚度的关系以及岩石圈有效弹性厚度各向异性等问题进行初步探讨。  相似文献   

水平螺旋槽管壁面液膜传热特性的研究   总被引：2，自引：0，他引：2  

武雅洁  梅宁  李艳 《中国海洋大学学报(自然科学版)》2005,35(1):167-172

对管式高效传热传质元件的工程应用进行理论探讨，根据水平螺旋槽管壁面液膜的形成机理及传热特性，通过建立单组分流体的物理和数学模型，得到了液膜速度和厚度解析解，并分析了在蒸发、冷凝时水平螺槽管的表面几何形状对液膜厚度分布的影响。结果表明：液膜的厚度主要受表面张力和槽道表面几何形状的影响，在槽道内弯处较厚，而在槽道起始处较薄，冷凝时与蒸发时相比液膜厚度更薄，液膜分布更均匀。相对于光滑直管，水平螺旋槽管壁面液膜具有更均匀的厚度分布，具有更好的传热传质性能。  相似文献   

Preliminary optical classification of lakes and coastal waters in Estonia and south Finland     

Anu Reinart  Antti Herlevi  Helgi Arst  Liis Sipelgas 《Journal of Sea Research》2003,49(4):357

A preliminary optical classification of lakes in Estonia and south Finland which can also be used for small bays of the Baltic Sea is elaborated. The classification is based on the optical properties of water (diffuse attenuation coefficient, diffuse reflectance) and parameters that are routinely monitored in water bodies (Secchi depth, concentration of chlorophyll-a, total suspended matter and yellow substance). The data complex used for our classification covers different types of water ecosystems (ranging from oligotrophic to hypertrophic) and the variability of water constituent concentrations in the ice-free period in Estonia and south Finland. Using cluster analysis, we found 5 optical classes of waters: clear (C), moderate (M), turbid (T), very turbid (V) and brown (B). There is satisfactory correspondence between class of water, shape of diffuse attenuation coefficient and diffuse reflectance spectra and trophic state of the lakes.  相似文献   

Full-range equation for wave boundary layer thickness     

Ahmad Sana  Hitoshi Tanaka 《Coastal Engineering》2007

Full-range equation covering all the flow regimes in a wave boundary layer is proposed for the boundary layer thickness. The results are compared with the available experimental data and good agreement has been found. In case of wave boundary layers, there are three definitions of boundary layer thickness in use. Therefore, the full-range equation is derived for three of the definitions. The findings of this study may be useful in calculating suspended sediment transport in coastal environments and studying wave–current combined motion.  相似文献